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USER GUIDE FOR STRMLN : 

A BOUNDARY-LAYER PROGRAM FOR CONTOURED •• 
WIND-TUNNEL LINER DESIGN 


By 

E. Clay Anderson 
DCW Industries, Inc. 


SUMMARY 

This report is a user guide for a 2-D boundary-layer computer 

o 

code which was developed in order to process data for an 
arbitrary number of streamlines. Provisions are included for 
the computer code to determine either (1) mass transfer rates 
necessary for an effective boundary-layer displacement of 
zero thickness or (2) the effective displacement thickness 
for a specified mass transfer-rate distribution. The computer 
code has been developed to be compatible with other computer 
codes which are being modified and/or developed at the NASA- 
Langley Research Center in order to design the three-dimen- 
sional, contoured, wind-tunnel liner to be used in transonic 
testing of a laminar flow control (LFC) system installed on a 
supercritical airfoil section. 

Appendices to this report present a brief discription of the 
liner design procedure, representative liner calculations, 
adaptive-wall design for a two-dimensional wind tunnel test, 
and other applications. 

INTRODUCTION 

'--Th’e~ 'purpose of this report is to provide a user guide for a 
2-D boundary-layer analysis computer code (STRMLN) which is 
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compatible with the string of computer codes being modified 
and/or developed at the NASA-Langley Research Center in order 
to design a contoured wind-tunnel liner. The program is 
written in FORTRAN IV programming language for use on the 
CDC CYBER series computers at NASA-LRC. A brief description 
of the current design procedure is given in Appendix A. 

Program STRMLN is a modification set for program LTBPG . 

LTBPG ' is * identical totthepprogramudiscussed in reference" 1 
except that LTBPG is restricted to the analysis of perfect 
gases having a constant or an effective-constant ratio of 
specific heats. The formulation of the boundary-layer 
equations and the numerical solution method are discussed 
in reference 2. In the present report, it is assumed that 
references 1 and 2 are available, and therefore, only the 
input data necessary for the present applications are dis- 
cussed. A number of irtput paramaters discussed in reference 
1 have been assigned values appropriate for two-dimensional 
flow calculations. One should see reference 1 for the defini- 
tion of paramaters not given in this report or if program 
STRMLN is to be applied to other classes of boundary-layer 
flows . 

Representative applications of program STRMLN are presented 
in appendices B, C, and D. A second modification set for 
program LTBPG which interacts the boundary-layer displacement 
with the mean inviscid flow field in nozzles is briefly dis- 
cussed in appendix E. 

PROGRAM DESCRIPTION 

The modification set (STRMLN) has been developed for program 
LTBPG to analyze the boundary layers along an arbitrary 


number of streamlines. The primary purpose of this mod- 
ification set is to provide the viscous-displacement or the 
effective-displacement (reference 3) corrections along - 
streamlines defining the contours of a non-porous wind-tunnel 
liner. The liner is to be used for transonic testing of a 
laminar flow control (LPC) system on a large-chord yawed- 
wing model. However, the resulting computer code can be 
applied without modification to the design of adaptive-wall 
two-dimensional wind tunnels or other applications where it 
is necessary to analyze a large number of streamlines. The 
analysis is restricted to flow fields where the curvature- 
induced crossflow-velocity components remain small. 

For the design of a contoured wind-tunnel liner, more than 
200 streamlines may be needed to adaquately describe the 
liner geometry. This analysis may have to be repeated in 
order to arrive at a configuration which is compatible with 
test requirements and constraints imposed by the existing 
wind-tunnel geometry. Because of the large amount of data 
to be processed, it is essential that the computer code be 
relatively free of data manipulation by the user. The present 
computer code has been applied to representative liner-design 
calculations and, with the exception of those streamlines 
where it is necessary to apply suction to avoid flow separa- 
tion, the data can be processed automatically . For the 
calculations which have been made to date, there are always 
a few streamlines at the airfoil/liner junction (both above 
and below the model) which must be treated individually. 

The procedure used to define the suction-rate distributions 
necessary to maintain attached boundary-layer flow requires 
the user to specify the distribution for each of these stream- 
lines and repeat the procedure until a satisfactory solution 
is obtained. In order to accelerate this cut-and-try process, 

. ...,„t.he ...c.o rapu-ter code contains logic to permit a number of 


3 


suction-rate distributions to be specified for a given stream- 
line. The distributions which result in attached boundary- 
layer flow are monitored to select the one for output which 
has the least total mass transferred from the channel. A 
numerical procedure is being developed to permit the computer 
code to specify the suction-rate distributions automatically; 
but, to '"date , this . procedure- -has' not been satisfactory for’ the 
abrupt adverse pressure gradients associated with the leading 
and trailing edges of the airfoil. This procedure provides 
satisfactory results for relatively gradual adverse pressure 
gradients but tends to concentrate high levels of suction 
over a short distance if the pressure gradients are abrupt. 

The procedure for locating the streamlines for which a flow 
separation is predicted is discussed in the following paragraph 

For the present application, the basic liner shape may be 
considered as rectangular, and data files are generated for 
each of the four surfaces. These files contain the necessary 
geometric and pressure coefficient data for approximately 
50 streamlines defining each liner surface. The boundary- 
layer solutions along the streamlines forming the liner walls 
above and below the airfoil surface can be obtained without 
difficulty. For the data files containing the streamline 
data defining the airfoil/liner junction (i.e»s, the endplates), 
those streamlines which require suction generally pass within 
a distance of +-6 percent chord above and below the junction 
line, respectively. To determine the number of streamlines 
which need to be treated individually, one first inspects the 
inviscid flow-field data and determines the streamline index 
associated with the streamline which is located approximately 
10 percent chord above the airfoil surface. Next a job is 
submitted to solve the streamlines in sequence up to the pre- 
viously determined index. Provisions are included within the 
computer code to monitor the solutions so as to terminate 
those along streamlines for which flow separation is predicted 
— ~ and, J a.d.vanoe-"'to*“t'h'e "next streamline. Printed messages identiy 



the streamlines for which it is necessary to provide suction- 
rate distributions in order to keep the boundary-layer attached. 
These streamlines are then computed individually (one for each 
job submitted) until all streamlines requiring suction have 
been satisfactorily solved.' The remainder of the problem 
can then be solved without difficulty. 

Provisions are also included in the program to permit the 
analysis of boundary layers on the airfoil surface with 
options to compute the effective airfoil shape or to compute 
mass-transfer distributions which result in an effective dis- 
placement thickness of zero. 

Representative applications of the computer code are presented 
in appendices B, C, and D. Appendix B summarizes the numerical 
results obtained for an adaptive two-dimensional wind-tunnel 
wall shape and presents comparisions with the experimentally 
determined shape. Applications of STRMLN to the analysis of 
boundary layers on the airfoil surface are presented in 
Appendix C. A representative liner calculation is discussed 
in Appendix D. 


INPUT DATA 

Input data to program STRMLN is via punched cards and disc 
file. The punched card data in namelists and file data are 
presented in the order which they are read by the program. 

The namelist name or disc file unit number is given first 
and is followed by a description of the data read. References 
to input data not previously defined will be found in a sub- 
sequent alphabetical listing. The description of data which 
is referred to and not/' defined in the present report is given 
in reference 1. Many of the input variables which are dis- 
cussed in reference 1 have been assigned values appropriate 
for the present application; only those variables which are 
chan ged frequently are described herein. 



NAMELIST/P IXPT/* 


IADW 


IBL1 


IBL2 


ICQZDST 


IE 


INLT1 


JJI 


JJM 


Indicator for thermodynamic condition at wall. 

= 0, non-adiabatic 
= 1, adiabatic 

Subscript in the SS array defining the point 
where the boundary-layer solution is initiated. 
Note that for INLT1 / 0, SS(IBLl) must have 
the same value as x read from file NT.4 . 

Subscript in the SS array defining the point 
where the boundary-layer solution is terminated. 
Program dimensions require 

5 < IBL2 - IBL1 + 1 <101. 

Indicator, for mass transfer at wall. If non- 
zero, the mass transfer distribution for an 
effective displacement thickness of zero is 
determined for the interval. 

SS ( IBL1 ) < x < SS (IBL2 ) 

Number of grid points in the surface normal co- 
ordinate used to solve the boundary-layer 
equations. Maximum permitted value is 101 and 
is the value normally used. 

Indicator. Identifies the record on file NT.4 
to be used as the boundary-layer starting 
solution when all streamlines to be solved 
have the same upstream environment. 

Index which identifies the streamline at which 
the boundary-layer solutions begin. A value of 
zero is set to 1. 

Indicator giving the total number of streamlines 
to be analyzed by the boundary-layer computer 
code. A value of zero is set to 
JJM = JIMAX - JJI + 1 

If IMCQD^O, JJM is assigned a value of 1. 
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KST2 


KTRANS 


KTRNSN 


LAMTRB 


Indicator used in conjunction with the output 
file NT2. A zero value indicates that no pre- 
vious records have been written on file NT2 
and JJI is assigned a value of 1 unless NT2 
is assigned a value of 0. If a value of 1 is 
assigned to KST2 and NT2 is non-zero, records 
on file NT2 are read up to JJI-1. The data on 
file NT2 must be recorded in'-- sequential order. — 

Indicator for boundary-layer transition. If 
non-zero, a continuous transition from laminar 
to turbulent flow is computed. LAMTRB is set 
to 1 and CHICRT must be defined. 

Subscript in XSTA array defining the point at 
which the boundary-layer is assumed to make the 
transition from laminar to turbulent flow in- 
stantaneously. A value of 0 is set to IBL2- 
IBL1 +1. If a non-zero value is input, LAMTRB 
is set to 1. 

Indicator for type of boundary layer. 

= 1 for laminar or transitional flows. 

= 2 for turbulent flow. 

Note that KTRANS or KTRNSN ? 0, must have 
LAMTRB = 1. 


NASY - Indicator for previous file processing by assembly 

program. 0 indicates that the input data file 
(NT1) has not been processed by the assembler 
program (see Appendix A). A value of 1 indicates 
input data has been processed by assembler 
program. 


NT1 
NT2 
NT 3 


Unit number for streamline data input file. 
Unit number for streamline data output file. 

t 

Unit number for auxiliary output file. 
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NT4 - Unit number for boundary-layer starting solution 

input data file. 

Note: NT1 must be defined. NT2, NT3, and NT4 are bypassed if 

assigned a value of zero. 


NAMELIST/FLTPT/ 

BO - If IADW is Oj the surface temperature is held 

constant at the value BO*TSTAG. 


CHICRT 


CHORDWN 


DX 

DXMAX 

G 

PRL 

PRT 

PSTAG 


Vorticity Reynolds number at which a continuous 
transition from laminar to turbulent flow is 
initiated. KTRANS must be non-zero for a con- 
tinuous transition, and KTRNSN is assigned a 
value of 0. The normal range for the critical 
vorticity Reynolds number at the transition 
point is 

2000 < CHICRT < 4000 

for subsonic 2-D flows. If KTRANS is non-zero 
and CHICRT is 0, transition is initiated at a 
momentum thickness Reynolds number of 1000. 

Chord length (in feet) measured in the wing- 
normal direction (i.e. normal to the leading 
edge of a yawed wing). 

Initial streamwise integration stepsize. 

Maximum streamwise integration stepsize permitted. 
Ratio of specific heats. 

Molecular Prandtl number. A value of 0 is set to 

0.7. 

Turbulent Prandtl number. A value of 0 is set 
to 0 . 9 • 

Stagnation pressure, PSIA. 


8 



RECHSTL - Reynolds number based on the streamwise chord 

length. The value is input as 0 and is com- 
puted by the computer program. 

RECHWN -- Reynolds number based on the wing-normal chord 

length. The value is input as 0 and is computed 
by the computer program. 

TSTAG - Stagnation temperature in degrees Rankine . 

NAMELIST/MASSTR/ 

ICQD - A non-zero value indicates that the mass transfer 

rate distribution changes discontinuously . Up 
to 10 discontinuous changes in the mass transfer 
rate can be entered in the array CQD. 

IMCQD - If non-zero, the streamline JJI is solved any 

desired number of times using different mass- 
transfer rate distributions. The present name- 
list is read after each solution is obtained, 
and the solutions are terminated when an end- 
of-file is read. 

IMT - If non-zero, a mass transfer rate distribution 

or distributions must be specified. 

KEMT - Array of subscripts for the XSTA array (see 

reference 1) indicating the points where mass 
transfer is terminated or changes discontinuously. 
KEMT is dimensioned for 10 values. A shifted 
index is used to define the XSTA array and is 
defined as follows: 

XSTA (N) = SS (K) 
where 

N = K - IBL1 + 1 
and 

IBL1 < K < IBL2 
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KIMT 


Array of subscripts in the XSTA array (see 
reference 1) defining the points where mass 
transfer is initiated or changes discont inuously 
from one level to another. KIMT is dimensioned 
for 10 discontinuous changes in the mass transfer 
rate. ■•-If ICQD- is .non^.zero, the array CQD must 
have the same number of entries as the arrays 
KIMT and KEMT. If the mass transfer rate dis- 
tribution has a continuous variation, only 1 
value of -KIMT . and KEMT are entered, and the 
mass transfer rate distribution CQS must be en- 
tered at all points in the XSTA array. Note 
that CQS can change from 0 to a finite value or 
from a finite value to 0 at the points where 
mass transfer is initiated or terminated. 

CQD - Array of mass transfer rates which change dis- 

continuously from one level to another at the 
beginning of intervals defined by KIMT and KEMT. 
ICQD must be non-zero and the same number of 
entries must be made for KIMT, KEMT, and CQD. 

CQD is dimensioned by 10. 

CQS - Array of mass transfer rates which have a con- 

tinuous variation in the interval KIMT (1) to 
KEMT (1). Outside this interval the values of 
CQS must be entered as zeros. The number of 
entries in the CQS array is IBL2-IBL1 + 1 <_ 101. 

DISC FILE INPUT DATA. UNIT NT1 . 

The disc file input data contains a number of variables which 
are not used in the boundary-layer computer code. These data 
have been used in previous processing steps or will be used in 
the post processing step. A description of all input data is 
given, and the -data which is necessary for the boundary-layer 
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calculations are superscripted with an asterisk- (?) . In the 1 

description of the disc file input data, the records and the 

i ‘/i 

variable names are presented in the order,/ ik- which they are ' read 
and the format is given. The variable names are then defined 
in alphabetical order following each record. 


DESC2 
DESC2 , 

DESC1 
DESC1 , 

DESC3 

DESC3, 

AMIN, 

ALAMD 

AMI* 

AMIN 

CPST 

QI 

QIN 

IMAX*, 

ILE* 


IMAX* 




(8A10) 

Alphanumeric description data identifying one 
previous processing procedure. 

(8A10) 


f. 

i 


Alphanumeric description data identifying one 

previous processing procedure. 

(8A10) (Read only if NASY = 1) 

Alphanumeric description data identifying 

previous processing procedure. 

QIN, ALAMD, AMI*, QI, CPST (8E16.8) 

- Sweep angle in degrees. 

- Freestream Mach number. 

- Component of freestream Mach number normal 
to the wing leading edge. 

- Pressure coefficient at sonic condition. 

- Freestream velocity, feet per second. 

- Freestream velocity component normal to the 
wing leading edge, feet per second. 

ILE*, ITE*, JIMAX*, JIL*,OTU*, (1615) 

- Subscript in the SS array identifying the 
location of wing leading edge. For boundary- 
layer calculations on the airfoil surface JJI 

is input as JIL or JIU; IBL1 = ILE-1; IBL2 ='ITE-1 
NT2 = 0; and JJM is either 1 or 2 . 

- Related to the number of data points in the 
SS array. IMAX is used by the programs gen- 
erating the streamline data and two of the data 
points are not recorded on the input file. 
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The boundary-layer computer code uses IMAX1 = 

IMAX-2 and is dimensioned for 5 IMAX1 <_ 101. 

ITE* - Subscript in the SS array identifying the 

location of the wing trailing edge. 

JIMAX* - Number of streamline records written on this 
disc file. 

JIL* - Index identifying the streamline associated 

with the airfoil lower surface-.- 
JIU* - Index identifying the streamline associated 

with the airfoil upper surface. Note ahead 
of and behind the airfoil, the y-values associ- 
ated with these two indices are the same. 

ALSBD, ALSTD , YJSB, YJST, SPAN, YSTLC (8E16.8) (Read 
only if NASY = 1) 

These quantities are not used by program STRMLN. These are 
parameters pertaining to the assembly of streamline data into 
either a "flat" end plate (wall) or an "octagonal" endplate 
(wall). They must be read and transferred through this program 
for use later in the processing step. 

JSB, J ST (1615) (Read only if NASY = 1) 

Above remarks apply to these quantities also. 

XS(N), YS (N) , ZS (N ) , CPS* (N ) , SS*(N), DSTREF* (N ) , .. - 'f 
r'l ' N =Ml s i IMAXl - (;8E16 &.) 

CPS* - Pressure coefficient on streamline. 

DSTREF* - Effective displacement thickness normalized 

with respect to the wing-normal chord length. 

Input as zero to initialize the array. 

SS* - Distance along the streamline normalized with 

respect to the wing normal chord length. 
r-„ 

XS,YS,&ZS - Carte.sian coordinates of streamlines defining -***•■- 

4 ‘-v 

the inviscid liner contours. The coordinates 
• - - r • ■ are normalized with respect to the wing-normal 

chord length. These coordinates are corrected 
for the displacement effect in the post pro- 
cessing step. 
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DISC FILE INPUT DATA. UNIT NT^ . 


Input data file NT4 is used if all streamlines to be solved 
have the same upstream environment. The boundary-layer 
starting solution data written on unit NT4 is generated using 
program LTBPG. 

INLT - (Unformatted) 

INLT - Record index used in conjunction with INTL1 

to select the desired boundary-layer starting 

solution . 

ff> ' 

X, XI, Z, Z0L j RO’, BETA, PP, NIT, K, TE, UE, XM, R0WE, XMUE, PE, 
DUEDS, CF, QW, HG, EPSVD, DM, KEP, NITT0T, CHIMAX, GAMMA, XIBAR , 
REX, QD0T, HG1 , HG2 , STE, STINF, CHEDGE, CH, CHREY, HAFCF , CH0CF , 
X0REFL, Z0REFL, R0REFL, DEL0X , TH0REF , DST0RF , DSAX0R , TH0DEL , 
DST0DL, DST0TH, DSTRAX , DSTARK , DELST, DEL, THET, RETHET, XN, 

YY, Y0VTHT, FC, FCN , F2NN, EPSPL, A0BP,"'tC, TCN, TH, R0R0E , C, 

CP, CHI, ETAINF, PNC, DS, XI0LD, X0LD, UER02, R0WEP, XMUEP, CFI, 
AIB, Y, FCP , THP , DN, FI, FIN, FINN, F2, F2N, Tl, T1W, T1NN, 

T2 , T2N, T2NN . - (Unformatted) 

Definitions of the variables appearing in this record are given 
in reference 1. 


OUTPUT DATA 


Output data from program STRMLN is via printed computer forms • 
and disc files. The printed output is primarily a listing of 
some of the input data, those boundary-layer results needed for 
determining mass-transfer distributions along individual stream- 
lines and error messages. The disc file output is the informa- 
tion which is transferred to the post processing step in order 
to make a wall correction for the effective displacement thick- 
ness . 




13 



PRINTED OUTPUT 


The printed output provides only a minimum amount of informa- 
tion. The previously defined input is listed first and is 
followed by a listing of the geometric coordinates and edge 
conditions for the first streamline to be solved. Geometric 
coordinates and edge conditions for subsequent streamlines are 
not listed. Solution data, for each streamline consist of 
a listing of the surface mass-transfer rate (if applied), dis- 
placement-thickness, and effective displacement-thickness (if 
surface mass transfer is applied) distributions. These data 
are followed by a message identifying the streamline. 

DISC FILE OUTPUT DATA. UNIT NT2 . 

Unless specified, all data written on the output file have been 
previously defined, and only the variable names and the format 
for each record is given. 

INF0 ( 1 ) , INF0 ( 23)3 INF0 ( 22 ) ( YSTLTB* , A7 , 2X, 2A10*, 

from program STRMLN* ) 

INF0 - Processing data retrieved from the computer 

operating system when the job is first sub- 
mitted. INF0 ( 1 ) is the job name assigned by 
the operating system, - INF0(23) is the date, 
and INF0(22) is the time of execution. This 
information is put on the file by every pro- 
gram which processes it. 

DESC2 - (8A10) 

DESC1 - (8A10) £ 

DESC3 - (8A10), written only if NASY = 1. 

AMIN, QIN, ALAMD, AMI, QI, CPST (8E16.8) 

IMAX, ILE, ITE, JIMAX, JIL, JIU (1615) 

ALSBD, ALSTD, YJSB, YJST, SPAN, YSTLC ( 8El6 . 8 ) (written 

only if NASY = 1) 
(written only i,f NASY = 1) 



/ 

BO', CH0RDWN, PSTAG, RECHWN, RECHSTL, TSTAG (8E16.8) 
IADW, IBL1 , IBL2 , JJI, JJM, KTRANS, KTRNSN, LAMTRB (1615) 

The 10 records listed above are written only when the job is 
first submitted. 

XS (N) , YS(N) , ZS (N ) , CPS(N), SS(N), DSTREF(N), N =1, IMAX1 
(8E16.8) 

AUXILIARY DISC PILE OUTPUT DATA. UNIT NT3 . 

The auxiliary output file is used for plotting purposes only, 
and the data written is variable. 
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APPENDIX A 


OUTLINE OF THREE-DIMENSIONAL WIND TUNNEL LINER DESIGN PROCEDURE 

FOR THE NASA LFC EXPERIMENT 

Currently available inviscid transonic and boundary-layer 
analysis computer codes are being modified at the NASA-Langley 
Research Center in order to design a non-porous wind-tunnel 
liner for the purpose of testing a laminar flow control (LFC) 
suction system installed on a super-critical airfoil section. 

The proposed test is to be conducted in the Langley 8-foot 
transonic pressure tunnel at a freestream Mach number of 0.82. 

The streamwise model-chord length of 7 feet results in a very 
small tunnel height-to-chord ratio. This experiment and test 
facility are being designed to simulate free-air flow about 
an infinite aspect-ratio yawed wing, so that two-dimensional 
design and analysis methods are applicable. For this case, the 
three-dimensional flow field is obtained by the addition of a 
constant crossflow-velocity component to the two-dimensional 
flow-field solution obtained for the airfoil at the appropriate 
reduced freestream Mach number. The streamlines (particle paths) 
forming the liner contours in the test section are to be 
determined from this velocity field. In order to make a smooth 
connection with the existing tunnel upstream of the contoured 
test section, the liner contour is obtained by superimposing the 
three-dimensional airfoil (yawed wing) flow-field perturbation 
onto the velocity field of an appropriate wind-tunnel contraction 
section. 

The sensitivity of high-speed channel flows to minor variations 
in the effective area-ratio distribution requires that viscous- 
displacemant corrections be made to the inviscid liner contours. 

A two-dimensional boundary-layer analysis is appropriate for all 
streamlines forming the liner contour except those near the 
airfoil/liner junction, since outside of this region, the stream- 
,,-v..^..^-’“*Ti'ne curvature induced crossflow velocity components along the 
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APPENDIX A - Continued 


liner walls are negligible. The flow in the immediate vicinity 
of the airfoil/liner junction is three-dimensional and is to 
be controlled by applying variable suction rates on both the 
liner wall and the airfoil surface. 

The departure from conventional test procedures results in 
part from the necessity to (1) use a relatively large-chord 
model in order to provide adequate space for the LFC suction 
system, (2) satisfy Reynolds number scaling requirements, and 
(3) conduct the test in a quiet, low-turbulence transonic wind 
tunnel. For conventional testing at super-critical speeds, a 
wind tunnel height-to-chord ratio of 3 to 4 is necessary (see 
reference 4) to avoid excessive tunnel induced interference 
effects. In principal, the contoured-liner concept permits 
the use of tunnel height-to.-chord ratios of any size. However, 
at super-critical flow conditions, the liner is restricted in 
its application to a specified test model at a specified test 
condition and permits only limited variations in test conditions 
about the design point. 

There are four basic tasks in the contoured liner design pro- 
cedure and each of these tasks is accomplished by several 
computer codes which pass information along on disc files. 

These tasks are: 

(1) inviscid yawed-wing test-section design; 

(2) inviscid 3-D contraction design; 

(3) viscous displacement correction; and 

(4) post processing of data. 

In the following paragraphs, a brief description of each task 
in the liner design procedure is presented. 

INVISCID YAWED-WING TEST-SECTION DESIGN 

The airfoil ordinates which produce a specified set of performance 
characteristics at the design Mach number are determined using 
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APPENDIX A - Continued 


the NYU fast-solver computer code developed by Bauer, Garabedian, 
Korn, and Jameson, reference 5- The flow-field solution de- 
termined by this analysis is given in a transformed coordinate 
system, and the liner design procedure requires the flow-field 
properties and geometric data to be expressed in a Cartesian 
coordinate system. This could be obtained by an inverse trans- 
formaion of the coordinate system used in reference 5 - A more 
convenient procedure, however, is to recompute the flow field 
about the previously determined airfoil configuration using the 
TRANDES computer code developed by Carlson, reference 6, which 
determines the flow-field solution in the required coordinate 
system. 

The yawed-wing streamlines are determined by integration of 
the 2-D velocity components (i.e., those lying in a plane 
perpendicular to the wing leading edge) and the constant 
orthogonal sweep-velocity component. This 2-D integration 
starts from an initial set of ordinates well upstream of the 
airfoil and ends far downstream from it. A streamline assembly 
program forms the liner by translating these space curves 
according to sweep theory. The final program in this step 
interpolates these assembled curves onto a grid which is fixed 
in the tunnel in order to define the ordinates for starting 
the upstream integration through the contraction section, draw 
pictures of the test section lines, and perform other ap- 
propriate data processing. 

INVISCID 3-D CONTRACTION DESIGN 

c V 

To define the co-ordinates of the wind tunnel liner surfaces in 
the contraction section, it is necessary to superimpose the 
perturbation velocity field due to the yawed-wirjg onto the 
vel9 < c^ti£,s & .-f:0.rwan-a-X'iaM-y^ symmetric contraction which far up- 
stream can be related to the existing tunnel. The streamline 
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analysis is made using the General Electric, computer code 
STC (reference 7)- To provide the necessary data for super- 
imposing velocities and then integrating upstream, the outer 
boundaries of the approximate contraction must extend outside 
the boundaries of the existing wind-tunnel facility. The 
streamtube which closely approximates the far upstream area 
of the existing tunnel must closely approximate the area of 
the three-dimensional liner at the match point in order to 
conserve mass. An iterative-- process is necessary to arrive — 
at an approximate configuration which satisfies the requirements 
of the liner design in the test section and which contains the 
approximate streamtube that can be fitted within the boundaries 
of the existing wind tunnel. 

The velocity field for this axi-symmetric nozzlg solution is 
interpolated onto a cartesian grid and then a velocity per- 
turbation due to flow about the yawed-wing model is super- 
imposed. This produces a 3-D velocity field on a cartesian 
grid in the contraction section region. A 3-D integration 
upstream through this velocity-field gives shapes which are 
taken to be the inviscid 3-D contraction section. At the far 
upstream end, where the Mach number, is low, these lines are 
faired back into the existing tunnel. 

VISCOUS DISPLACEMENT CORRECTIONS 

The coordinates defining the physical liner surfaces are the 
inviscid liner coordinates plus the viscous displacement cor- 
rection or the effective displacement correction (reference 3) 
if surface mass transfer is applied along streamlines. The 
displacement corrections to the inviscid streamline coordinates 
and the suction rates necessary to maintain an attached boundary 
layer flow are determined using program STRMLN discussed in 
this report. If the displaced geometric coordinates of the 


19 



APPENDIX A - Concluded 


liner surface fall outside the existing wind-tunnel boundaries, 
it is necessary to repeat, some aspects of the two design tasks 
outlined above. Since the lines are faired into the existing 
tunnel far upstream, an approximate starting solution for the 
boundary layer is required. It is obtained from LTBPG. 

POST PROCESSING OF DATA 

The final step in the liner design procedure is to process the 
data to put it in a format suitable for engineering and 
fabrication purposes. The liner shapes are given as parametric 
space curves and need to be interpolated onto a suitable re- 
ference system showing elevation, cross-sections, etc. 

It can be seen that tasks 1 and 3 above are applicable to the 
design of adaptive wall two-dimensional wind tunnels. Results 
for such an application and comparison with the experimentally 
determined shape are given in Appendix B. 
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T" 

APPLICATION- SE STRMLN CODE IN 2-D 
ADAPTIVE-WALL WIND-TUNNEL DESIGN 


The numerical procedure being developed for contoured wind- 
tunnel liner design (see appendix A) has been applied to de- 
fine the wall shape for a two-dimensional adaptive-wall wind- 
tunnel test. The numerically determined wall shape is compared 
with the wall shape determined by Barnwell and Everhart* using 
an iterative analytical/experimental method. Complete details 
of this application are presented in reference 8, and only a 
summary of the results are presented here. 


The experimental data were obtained in the NASA Langley 6- by 
19-Inch Transonic Tunnel (reference 9). For this experiment, 
the slotted wind-tunnel walls were removed and replaced by 
nonporous flexible jack-supported plates as shown in figure 
B-l. The model used is an NACA-0012 symmetric airfoil section 
at zero incidence and the test conditions were: 


Freestream Mach number, = O .765 

Freestream Reynolds number, N R ro = 2.00x10 /cm 

3 |T 

Stagnation pressure, P Q ! = 1.42x10 Pa 
Stagnation temperature, T^. = 282 °K 



The chord length, c, is 15-24 cm (6 inches) and the model was 
centered at the (+) mark shown in figure B-l. Coordinates.- x/c 
and y/c . are measured along and perpendicular to the wind 
tunnel's horizontal plane of symmetry with origins at the cen- 
ter mark and the plane of symmetry, respectively. 


*A detailed description of the streamline procedure and results 
obtained from it have not yet been published by R. W. Barnwell 
and J.L. Everhart "of NASA Langley Research Center. The present 
analytical results are compared with their experimentally de- 
termined wall shape; we appreciate and acknowledge this early 
release of their data. In this paper their procedure and re- 
sults will be identified as Barnwell-Everhart . 
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NASA LANGLEY 6- BY 19-INCH TRANSONIC TUNNEL 
WITH SLOTTED WALLS REPLACED BY FLEXIBLE WALLS 


LT) 
I 

< 


O 

o 

ai 


< 

i 

< 

E ] o 

i — 

o 

LU 

oo 




<C LU 

O Cl 

o o 



t/1 LO 
LU LU 

O O 

z z 

LU LU 
Q Q 

D> © 


22 


Figure B-l. Schematic of two-dimensional "streamlined tunnel" experiment of Barnwell and Everhart. 
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The free-air solution about the airfoil geometry plus the viscous 
displacement correction were determined using the transonic 
design/analysis computer code discussed in reference 5 and the 
transonic analysis computer code presented in reference 6. At 
the upstream junction of the flexible and rigid tunnel walls 
(x/c = - 4.833) s 48 streamlines lying in the region 0 < y/c 
< I .583 were extracted from the free-air solution for boundary- 
layer calculations. The starting boundary-layer solution was 
taken as that for a two-dimensional model of the wind-tunnel 
geometry up to the junction of the flexible and rigid tunnel 
walls. From this starting solution, boundary-layer calculations 
were made along each of the streamlines to define the viscous- 
displacement thickness distributions. The two-dimensional 
boundary-layer equations cannot be integrated into the strong 
adverse pressure gradients associated with the streamlines 
near the airfoil/tunnel junction; it was necessary to make 
modifications to the pressure distributions in the vicinity of 
the airfoil leading edge. The procedure used in the analysis 
was to replace four values of the pressure coefficient along 
each of four streamlines lying within the region - 0.57 £ x/c — -4- 

f 

£ - 0.49 by the corresponding values associated with the first 
streamline for which the boundary-layer analysis did not predict — 
a flow separation. This in effect reduces the pressure gradients 
to a value slightly below that for a predicted boundary-layer 
separation. More drastic changes in the pressure coefficients 
than those used in the analysis had a negligible effect upon the 
solutions downstream of the leading edge. However, this procedure 
should not be used if the zone of the separation extends to an 
appreciable part of the sidewall. All of the modified stream- 
lines are close to the airfoil/tunnel-sidewall junction where 
the flow is three-dimensional and the present analysis is not 
expected to be very accurate either with or without the mod- 
ifications in this region. 
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Figure B-2 presents boundary-layer displacement-thickness 
distributions across the tunnel sidewall (or across streamlines) 
at six locations along the streamwise direction. In the region 
- 2.10 <_ x/c £ - 0.57j the distributions show the expected 
upstream influence of the model pressure field. The plot at 
x/c = - 0.57 shows the influence of a favorable pressure 
gradient upon the sidewall streamline far from the model. 

That is, the displacement thickness along these outer stream- 
lines at this location are less than those at an upstream 
location. The inflection point in the displacement-thickness 
distribution at x/c = - 0.31 is the result of different flow 
acceleration rates along the streamlines and a difference in 
distance measured along the. streamlines relative to their res- 
pective upstream pressure peaks. The termination point of this 
plot above the y/c-origin represents the airfoil thickness at 
this location. 

Displacement thickness distributions at x/c = 0.506 and 1.07 
show the influence of the trailing-edge compression and a 
representative downstream profile. The inflection point in 
the displacement thickness distributions at x/c = 1.07 is the 
result of the more rapid flow acceleration along streamlines 
near the corner of the (effective inviscid) blunt trailing 
edge of the airfoil. 

As shown in figure B-l, only the two lateral walls opposite 
the airfoil surfaces are flexible. Thus, at each streamwise 
location, x/c, the displacement thickness, 6*/c (x/c, y/c ) , 
must be integrated around the local cross-sectional perimeter 
and then be applied as a displacement of only the two lateral 
flexible walls. Since no displacement correction was made at 
the (upstream) junction of the fixed rigid wall and the flexible 
plates, the wall displacement correction is made relative 
(differential) to that at this junction. The total wall dis- 
placement is composed of two components: (1) the departure of 
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Figure B-2a. Calculated boundary layer displacement thickness distributions on the tunnel sidewall 
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Figure B-2b." Concluded. 
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the outer bounding streamline from a straight line - the com- 
pressible blockage correction and (2) the boundary-layer dis- 
placement on all of the wind-tunnel walls - the viscous block- 
age correction. Comparisons of the numerically determined wall 
shapes with the experimental data of Barnwell and Everhart are 
shown in figure B-3- The symbols show the total experimentally 
determined displacement at the 11 jack-point locations indicat- 
ed in figure B-l. The maximum difference between the analytical- 
ly-determined wall-displacement corrections and the experimental 
data is at x/c = - O. 667 . At this location, the displacement 
corrections differ by less than 8 percent and the ratio of the 
two coordinates, which indicates the difference in the tunnel 
area ratio, is 1.002; The displacement contributed by the 
deflection of the outer bounding free-air streamline is also 
shown in the figure. 
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Figure B-3.“ Comparison of flexible nonporous wall shapes for an NACA 0012 
airfoil test at zero lift in the Langley 6- by 19-inch transonic tunnel. 
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APPENDIX C 


APPLICATION OF STRMLN CODE TO ANALYSIS OF 
SUCTION ON AN LFC AIRFOIL 


Surface mass-transfer rate distributions required 'to maintain 
laminar flow on an LFC airfoil designated as YNRB-12-2-77 
(Reference 10) have been determined using a laminar-flow 
stability analysis. The present two-dimensional computer code 
has been used to determine the effective displacement-thickness 
distributions on the airfoil surface. The effective displace- 
ment thickness, A#/c, is defined as (see Reference 3): 

f X 

A*/c = 6 */c + — ----- CQ(x/c ) d(x/c) (1) 

e e J 0 


where 



oo 

r- -, 

S*/c = 

• 


pu 

0 . 

1 P U 

e e 


d(y/c) (2) 


In these expressions, pu and p g u e are the tangential mass fluxe 
at a local point and at the boundary-layer edge, respectively; 
and CQ(x/c) is the surface mass-transfer rate distribution 
normalized with respect to the freestream mass flux. (See 
Reference 1.) 

The assumed test conditions for this LFC airfoil calculation 
are as follows : 

Freestream Mach number, M = 0.891 

00 6 

Freestream Reynolds number, N^ g m = 15x10 , based on the 

streamwise chord 

Stagnation pressure, P = 5-87 PSIA 
Stagnation temperature, Tc = 508° R 
Wing-Normal chord, c = 6.553 ft. 
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Streamwise chord, c = 8.0 ft. 

. Sweep" angle, A = 35° 

The airfoil ordinates and surface pressure distributions are 
shown in figure C-l and the surface mass transfer rates are 
shown in figure C-2. Figure C-3 presents the corresponding 
effective displacement thickness distributions for laminar 
flow. The effective displacement thickness is negative 
over approximately 80 percent of the chord on the lower sur- 
face and over approximately 25 percent of the chord on the 
upper surface. The effective displacements for this case are 
not large enough to result in a significant change in the air- 
foil pressure distribution; however, in order to check the steps 

w '• • 1, 

in the design procedure, the displacement thickness was^ added 
to the airfoil ordinates, and the free-air flow field about the 
effective airfoil geometry was recomputed using the transonic 
analysis computer code discussed in reference 5 • Figure C-4 
presents a comparison of the pressure distributions and sonic 
line locations for the airfoil ordinates and the effective 
airfoil ordinates. 

On the outboard portions of the airfoil model, the turbulent 
liner-wall boundary layer will contaminate the flow over the 
airfoil surface and result in turbulent zones of flow. To 
determine the effective turbulent displacement-thickness dis- 
tributions in flow regions where the laminar-stability mass- 
transfer rates are applied, the boundary layer was assumed to 
be fully turbulent, and the resulting effective displacement- 
thickness distributions are shown in figure C-5. These tur- 
bulent A*/c's are sufficiently large to significantly alter the 
airfoil pressure distribution. This influence can be controlled 
by subtracting the effective displacements from the airfoil or- 
dinates but would result in a compound airfoil surface. An 



□ - Upper surface 
O- Lower surface 



Figure C-l.- 2-D pressure coefficient distributions and 
coordinates for LFC Airfoil YNRB-12-2-77. 







Figure C-2.- Surface mass-transf er rate distributions for 

laminar flow stability on LFC Airfoil YNRB-12-2-77 
at test conditions. 
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Figure C-5.- Effective displacement-thickness distributions 
for turbulent flow with laminar flow stability 
suction rates applied on LFC Airfoil YNRB-12-2-77 
at test conditions. 
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alternate procedure is to determine mass-transfer rate dis- 
tributions for an effective displacement of zero by an iterative 
solution of equation (1) and use this solution as a guide to 
specify mass-transfer rates which result in displacements having 
a negligible influence on the airfoil performance. The mass- 
transfer rate distributions for zero A#/c are presented in 
figure C-6. For these calculations, convergence was assumed 
to be |A*/c|<_ 10 . The rapid growth in CQ near the leading 

edge is the result of assuming fully developed turbulent flow, 
and the downstream irregularities reflect changes in the pres- 
sure gradient. (See figure C-l.) By inspection of figure C-6, 
a relative mean constant CQ of - 1.5x10 was selected and ap- 
plied to both the upper and lower surfaces of the airfoil. The 
resulting effective displacement thickness distributions are 
shown in figure C-7. For this airfoil at these test conditions, 
the constant value of CQ appears to be satisfactory. 
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Figure C-6.- Calculated surface mass-transfer rate distributions 
to maintain zero effective displacement thickness 
for turbulent flow on LFC Airfoil YNRB-12-2-77 at 
test conditions. 
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Figure C-7.- Effective displacement-thickness distributions 

for turbulent flow with a constant suction rate, 

CQ = - 0.0015, applied on LFC Airfoil YNRB-12-2-77 
at test conditions. 
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REPRESENTATIVE LINER CALCULATIONS 

Some of the computer codes needed for a complete inviscid liner 
design are currently being developed by NASA. However, the 
boundary-layer computer code developed for this application 
has been applied to the streamlines forming the contours of a 
liner test section for a preliminary LPC airfoil designated as 
971T. The velocity-field superpositioning used to obtain the 
liner contraction, section contours (see appendix A) has been 
bypassed; streamlines and the local flow-field properties along 
them have been determined for free-air flow about the airfoil 
at the following conditions: 

Freestream Mach number, M„ = 0.82 

00 6 

Freestream Reynolds number, N Rg = 26x10 , based on the 

streamwise chord 

Stagnation pressure, Po = 13 PSIA «. — l-~ 

O Q ^ » 

1 Stagnation temperature, Tp = 508 R 

o' 

Wing-normal chord, c = 6.49 ft. 

Streamwise chord c = 7-0 ft. 

Sweep angle, A = 22° 

The determination of suction-rate distributions required to 
maintain an attached boundary layer on the liner sidewalls near 
the airfoil/liner junction is the only part of the design pro- 
cedure which is not fully automated. For this part of the an- 
alysis, it appears that the most satisfactory method is a cut- 
and-try approach. That is, one specifies a number of suction- 
rate distributions along each streamline for which the flow 
separates in the absence of suction. The computer program 
monitors the solutions corresponding to suction-rate distributions 
which result in attached boundary-layer flow, and selects for 
output the one for which the total mass removed is least. For 
the present calculations, the suction-rate levels have been re- 
stricted to'' the range"- 0.008 £ CQ £ 0 in increments (multiples) 
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of 0.002; a maximum of three different levels of CQ was con- 
sidered for a given distribution. However, as many as 10 
discontinuous changes in the suction-rate level may be specified 
in a given distribution, or the distributions may be specified 
as having a continuous variation. 

The suction rates indicated in figure D-l were obtained by 
specifying approximately 10 different suction-rate distributions 
for each of the six streamlines that required suction to maintain 
attached flow. The figure shows the ordinates of the airfoil 
and the suction-rate distribution on the liner sidewall near 
one of the airfoil/liner junctions. The suction rates ahead of 
the airfoil leading edge are shown to demonstrate that the com- 
puter code can be applied successfully to obtain attached boundary- 
layer solutions in the leading-edge region where the approaching 
boundary layer is relatively thick. In the LPC experiment, the 
boundary layer in the leading-edge region will be controlled by 
the use of much larger suction rates close to the leading edge. 

It is noted that the flow is three-dimensional in the junction 
region downstream of the leading edge where suction is applied. 

The suction-rate distributions specified by a two-dimensional 
analysis can be considered only as an indication of the rates 
necessary, and provisions will be included in the suction con- 
trol system of the experimental facility to vary these rates by 
a factor of two or more. 

A representative liner cross section taken downstream of the air- 
foil trailing edge is shown in figure D-2. For this liner the 
entrance cross section is rectangular. The steps in the liner 
walls result from different spanwise deflections of the stream- 
line that divides and passes through the different flow-field 
environments above and below the upper and lower surfaces of 
the swept-wing panel at lift. The streamlines forming the sur- 

.' 5 '”“- ' -»«• 

‘“■fa'ces of these steps originate at different locations along the 



- CQ = - 0.008 


V77m - CQ = - 0.004 
t 1 - CQ = - 0.002 



Figure D-l.- Suction rate distributions for attached turbulent 

boundary layer on contoured liner wall for flow about 
yawed-wing model. 


41 




Figure D-2.- Contoured liner cross section downstream from trailing 
edge of yawed-wing model. 
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span of the model at the attachment line instead of at an up- 
stream location on the wall liner. The boundary-layer dis- 
placement data are shown only for the streamlines that originate 
at an upstream location on the wall liner. For the cross section 
shown in figure D-2, the boundary-layer displacement correction 
is approximately 3*5 percent of the channel area and is represent 
ative of the displacement correction throughout the liner test 
section. For the high subsonic Mach number considered, the 
displacement is sufficient to result in choking unless the 
displacement corrections are made. 

Figure D-3 shows the boundary-layer displacement distributions 
near the airfoil/liner junction and the development of the step 
on the left liner sidewall at x/c = 0, 0.4, 0.8, and 1.2. Suc- 
tion causes a reduction in the boundary-layer displacement 
thickness in the immediate vicinity of the airfoil/liner junction 
and, at some streamwise locations, the resulting effective 
displacement corrections are negative. 
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INVISCID-VISCOUS INTERACTIONS IN WIND TUNNELS 

A modification set for the basic boundary-layer computer code 
LTBPG has been developed for interacting the viscous-displace- 
ment correction with the mean inviscid flow field in closed 
channels. This modification set (program CHANL) was developed 
in order to analyze the effects of a reduced test section area 
and/or suction upon subsonic diffuser performance. If there 
are no shockwave-boundary-layer interactions, the analysis is 
applicable to supersonic nozzles and wind tunnels. 

This program has been applied to analyze the diffuser performance 
of the NASA Langley 8-foot transonic wind-tunnel facility with 
a reduced test section area. A reduced test-section area in- 
creases the adverse pressure gradients in the diffuser and 
boundary-layer control may be required to maintain attached 
flow. The program contains provisions for analyzing the in- 
fluence of area distributed surface mass transfer upon subsonic 
diffuser performance. Interacted solutions are obtained by 
successive approximations using a weighted viscous-displacement 
geometry and area-ratio relations along with perfect-gas isen- 
tropic relations to define the inviscid edge conditions. The 
iteration procedure is continued until the maximum local change 
in the Mach number distribution between successive iterations 
is less than 0.1 percent or until a flow separation is predicted 
if a subsonic diffuser is being analyzed. For the analysis of 
subsonic wind-tunnel diffusers, the inviscid edge conditions 
corresponding to the diffuser's geometric area-ratio distribu- 
tion generally results in a predicted flow separation if the. 
test section Mach number is above about 0.5 and the diffuser 
half-angle is greater than about 2.0 degrees. For these cases, 
an inverse approach to the converged solution or a predicted 
flow separation is- utilized. This procedure requires that the 
first solution pass be determined using an approximate geometry 
'"f o r'~ WhTc h""'f 1 o w separation is not predicted. 
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Figures E-l and E-2' present representative results obtained 

for the Langley 8 -foot transonic wind tunnel operating at a 

Mach number of 0.80 without a test model in place. The ^ , ,-- s & 

stagnation pressure and temperature were 13 PSIA and -5 08°R . 

respectively, and the freestream Reynolds number was 4.0x10^/ 

foot for the calculations and two runs of the tunnel. The 

l 

experimental facility was operated at these conditions primarily 
to determine the noise level in the test section; velocity pro- 
file data were measured only at the first diffuser exit. Figure 
E-l shows the wind-tunnel geometry from the entrance of the test 
section to the end of the diffuser and the corresponding com- 
puted displacement geometries for two conditions: (a) no mass 

transfer and (b) a constant suction rate of - 0.001 applied 
. over the interval 56 <_ Z <_ 63 . The boundary-layer calculations 
were started at the entrance of the wind-tunnel contraction 
section and transition from laminar to turbulent flow was init- 
iated at the location where the vorticity Reynolds number first 
exceeded 2800. The behavior of the displacement geometries is 
as expected both with and without suction. The velocity profile 
computed for the case of no mass transfer is compared with the 
two experimentally determined profiles at the diffuser exit, 

Z = 131, in figure E-2. The experimental velocity profile 
data* were obtained using a pressure scanning probe. A time 
interval of approximately 20 seconds was required to scan all 
of the pressure probes used to determine a velocity profile. 

The agreement between the computed and experimentally determin- 
ed velocity profiles is only fair. However, the computed 
boundary-layer thickness and the displacement thickness differ 
from the experimental values by less than 10 percent. The 
same test conditions and tunnel geometry were analyzed using 
the computer code discussed in reference 7 and a flow separa- 
^^tip.n..j^as ! predict e<T upstream of the step in the diffuser. 


*The wind-tunnel experiments were conducted by Joseph Brooks 
of NASA'Langley'Researc.h 'Center , arid • the data' have not been 
published.'*' The release of these data is appreciated and 
acknowledged . 
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APPENDIX E— Concluded 

The two solution procedures showed excellent agreement with 
each other up to the entrance of the diffuser. Additional 
comparisons with more complete experimental data taken at 
several locations along a diffuser axis are necessary to 
establish the accuracy of the numerical solution procedures. 

Provisions are included in this modification set to permit . 
the analysis of tangential slot injection, and this option has 
been debugged using an over-simplified turbulence model. A 
more complete turbulence model has not been included since 
this method for controlling the boundary layer in the diffuser 
was considered to be incompatible with the low noise level 
requirements for laminar flow control experiments. 
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